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We study a large- TV generalization of J1-J2 Heisenberg model on square lattice - an Sp{2N) spin 
model. The possible quantum spin liquid phases of the Sp{2N) model are studied using the SU (2) 
projective construction. We find several spin liquid states at least in large A'' limit, which include 
SU{2) TT-flux state, 5*17(2) chiral spin state and Z2 spin liquid states. All those spin liquid states 
have non-trivial quantum orders. We show how projective symmetry group, which characterizes 
quantum order, protects the stability of even gapless spin liquids. We also study the continuous 
quantum phase transition from the SU{2) 7r-flux state to the SU{2) chiral spin liquid state, and 
from the SU (2) tt-Aux state to the Z2 spin liquid state. We show that those phase transitions are 
beyond the paradigm of Landau symmetry-breaking theory. The first phase transition, although a 
Z2 symmetry breaking transition, contains critical exponents that are different from those obtained 
from the Ginzburg-Landau theory of Ising universality class. The second transition does not even 
involve a change of symmetry and has no symmetry breaking order parameter. 
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I. INTRODUCTION 

The SU{2) projective construction (or, more generally, 
slave-boson theory [1, 2]) for 2D spin liquids in spin-1/2 
systems has been studied over ten years [3, 4]. Upon dop- 
ing, the SU (2) projective construction provides a quite 
complete theory for underdoped high Tc superconductors 
[5, 6]. The theory explains the strange Fermi surfaces 
[5, 6], the strong (tt, tt) spin fluctuations [7, 8], and the 
temperature dependence of superfluid density [9, 10] for 
underdoped samples. 

At half filling, the SU{2) projective construction can 
also be used to describe various spin liquids. In fact, 
the SU (2) projective construction (or more generally, the 
slave boson theory) have predict many different spin liq- 
uids, such as the algebraic spin hquid [8, 11, 12], chiral 
spin liquid [13, 14], Z2 spin liquids [15-18]. It was also 
shown that the SU (2) projective construction is capable 
of describing hundreds of different spin liquids that have 
the same symmetry but different quantum orders [12]. 

The above predictions of spin liquids and the classi- 
fication of spin liquids were based on mean-field calcu- 
lations. At the mean-field level, it is not very hard to 
design a spin Hamiltonian that realizes a spin liquid that 
has one of a few hundreds quantum orders. It is also not 
hard to find the mean-field ground state for a given spin 
Hamiltonian. The real issue is whether we should trust 
the mean-field results. It was argued [12, 16] that, if 
the obtained mean-field ground state is unstable (i.e. if 
the mean-field fiuctuations cause diverging interactions 
at low energies), then the mean- field result cannot be 
trusted and the mean-field state does not correspond to 
any real physical spin state. It was also argued that. 
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if the mean-field ground state is stable (if the mean-field 
fluctuations cause vanishing interactions at low energies), 
then the mean-fleld result can be trusted and the mean- 
fleld state does correspond to a real physical spin liquid 
state. 

However, the above statement about stable mean-field 
states is too optimistic. A 'stable mean-field state' does 
not have diverging fiuctuations at low energies. So it 
does not have to be unstable. On the other hand, it 
does not have to be stable either. This is because short- 
distance fluctuations, if strong enough, can also cause 
phase transitions and instabilities. Therefore, in order 
for a mean-fleld result to be reliable, the mean-field state 
must be stable {i.e. no infrared divergence) and the 
short-distance fluctuations must be weak. As we do not 
have any small parameters in the SU{2) slave-boson the- 
ory for spin-1/2 systems, the short-distance fiuctuations 
are not weak, even for stable mean-field states. Because 
of this, it is not clear if the mean-field results, even for 
the stable states, can be applied to the spin-1/2 model 
or not. 

In this paper, we will generalize the spin-1/2 model 
to a large- model. The large- A^ model can be solved 
approximately using the SU{2) projective construction. 
We will show that the SU{2) projective construction 
for the large- A^ model have weak short-distance fluctua- 
tions. Thus, the stable SU{2) mean-field states for the 
large- A^ model do correspond to real physical spin liquid 
states. The SU{2) mean-field results for the stable states, 
such as fractionalization, emergent gauge structures and 
emergent Fermi statistics, can be applied to the large- A^ 
model. 

We concentrate on a model with nearest-neighbor Ji 
coupling and next-nearest-neighbor J2 coupling. We find 
several mean-field phases that include SU (2) 7r-flux state, 
SU{2) chiral spin state and Z2 spin hquid state. 

The SU{2) TT-flux state is described by a low energy 
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effective theory that includes gaplcss Dirac fcrmions cou- 
pled to SU{2) gauge fields. Due to the non- vanishing 
interaction between the fermions and the gauge bosons 
down to zero energy, there is no free fermioninc or bosonic 
low energy quasiparticles in the SU{2) tt-Aux state. De- 
spite this, we show that the SU{2) 7r-flux state is a stable 
spin liquid state. It is a realization of algebraic spin liq- 
uids. [8, 12] 

The SU{2) chiral spin state and Z2 spin liquid state 
are both gapped and, thus, naturally stable. Both state 
carry non-trivial topological orders. 

Within our J1-J2 model, the quantum transition (i.e. 

the zcro-tcmpcrature transition) between the SU(2) tt- 
flux state and the SU{2) chiral spin state turns out to 
be a continuous transition. The transition breaks the 
time-reversal symmetry and has a well defined Z2 order 
parameter. However, we show that the critical proper- 
tics of the transition are not described by the 3D Ising 
universality class of the Ginzburg-Landau theory. This is 
because the transition not only break the Z2 symmetry, it 
also changes the quantum/topological order. The contin- 
uous quantum transition between the SU (2) 7r-flux state 
and the SU (2) chiral spin state is a new class of quantum 
transition. 

It has been shown that continuous quantum transitions 
are possible between two states with the same symmetry 
(but different topological orders). [19-23] The continu- 
ous quantum transitions are also possible between two 
states with the incompatible]^!] symmetries. [24]. In this 
paper, wc show that even symmetry-breaking transitions 
with well defined order parameters, sometimes are not 
described by Landau's symmetry breaking theory. So it 
appear that most quantum continuous transitions are not 
described by Ginzburg-Landau theory, regardless if they 
have symmetry breaking and order parameter or not. 

Wc also study the transitions from the SU (2) tt-Aux 
state to a Z2-linear spin liquid state, and from the SU (2) 
TT-flux state to a f7(l)-linear spin liquid state. The 
transtion between the SU{2) 7r-flux state to the Z2-lincar 
state is found to be a continuous quantum transition that 
does not break any symmetry. The transtion between the 
SU{2) TT-flux state to the ?7(l)-linear state, on the other 
hand; is a continuous quantum transition that breaks lat- 
tice rotaion and translation symmetry. What is surpris- 
ing is that the two seemingly very different transitions 
are described by the same critical point with the same 
set of critical exponents. 



II. FORMULATION OF SU{2) PROJECTIVE 
CONSTRUCTION 



First we would like to briefly review the SU{2) pro- 
jective construction. We will mainly follow the notation 
of Ref. [12, 25]. Let us start with spin-1/2 Heisenberg 



Model, and introduce the fermion representation of spins: 
i? = ^ Jij Si • Sj (1) 

Si = -^floc'^alifif} (2) 

If we regard H as an operator acting on the fermion 
Hilbert space, then we have enlarged the Hilbcrt space. 
Therefore we have to add extra constraints to reduce the 
enlarged Hilbert space to the original one for the spin 
system. The constraints are the one-particle per site con- 
straints: 

fiafia = 1 (3) 
^ fiafi/Sea/j = 0, fljlpCafi = (4) 

The two extra constraints in Eq. (4) arc results of the first 
one in Eq. (3). However in path integral formalism, if 
we enforce all constraints simultaneously by introducing 
some Lagrangian multipliers, a lattice gauge theory with 
SU (2) gauge group can be derived. Let us introduce 
some notations: 



(5) 
(6) 



These arc the only singlet bilinear forms of the pairing 
between site i and j. After some rearrangements, one 
has: 

Ignoring the irrelevant constant, the path integral La- 
grangian of the Heisenberg model Eq. (1) turn out to 
be: 



i 



1 + ^ 

2«0,i'?i' 



^"o,i(xli 



1) 



(ij) 



Vij 



xljXij) 



The constraint of one particle per site has been encoded 
by the Lagrangian multipliers in the second line: 
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If we do a particlc-holc transformation of the spin- 
down fcrmions together with a Hubbard-Stratonovich 
transformation, the Lagrangian can be written in a form 
with the explicit SU (2) gauge invariance: 



{ij) 



'-TrU 



(7) 
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where 



Ui. 



Xij Vij 
Vij Xij ■ 



(8) 
(9) 



The path integral that describes the spin- 1/2 system is 
given by 

Z = J D{,p)D{ayD{Uij)e'f''*'^ 



The SU (2) gauge transformation is given by 



i^i ^ V4 



^3 



W'- 



WiUijW] 



where all Wi e SU{2). 

There are some other useful relations. For each site, 
the conjugate of fundamental representation of SU (2) is 
equivalent to the fundamental representation itself. One 
therefore can introduce another SU{2) doublet (neglect- 
ing the site label): 



'tp = ia2il>* 



fl 



(10) 



Thus there are only three SU (2) gauge invariant bilinear 
forms for ^ fields on the same site 



s+ = 


= 




(11) 


s- = 




/i/T 


(12) 


s' = 












- flh) 


(13) 



They turned out to be the generators of the spin rotation 
symmetry. 

In the zeroth order approximation, or at the mean-field 
level, one assumes that the boson fields ■ and Uij get 
condensed, or more specifically, can be replaced by some 
time-independent c- numbers. In this approximation, the 
system is described by the following mean-field Hamilto- 
nian 



i 



{ij) 



(14) 



III. A LARGE-iV LIMIT OF SU{2) PROJECTIVE 
CONSTRUCTION 

The above mean-field approximation is a good one only 
when there are some reasons to suppress the fluctuations 
of the boson fields ^ and Uij. One way to suppress 
these fluctuations is to go to a large- limit. Actually 
the mean- field result is exact when N = oo. In this 
section, we try to answer the following questions: 

• What is the lattice spin model that the large- A/' 
limit corresponds to? 

• What is ground state of this lattice spin model? 

Here we present our answers of these problems first. 
Later we will see the detailed derivation of these answers: 

• The lattice spin model is a Sp{2N) spin model de- 
scribed by Eq. 18. The Sp{2N) spin model is a gen- 
eralization of the usual SU{2) spin model (which 
corresponds to A'' = 1 case). 

• The ground states of the model are Sp{2N) sin- 
glets. For small J2 (the next-ncarcst-neighbor cou- 
pling). The system is in an algebraic spin liquid 
state. [8, 12] The excitations are the fractionalized 
particles (spinous) coupled to SU{2) gauge fields. 
The gauge field is deconfined. For larger J2/J1, 
there is a continuous quantum phase transition and 
the system goes into the chiral spin state which 
breaks the time reversal symmetry. [14] 

It was debated for long time the existence of featureless 
Mott insulator if the unit cell does not have c;ven number 
of electrons. In that case, it seems that in order to have a 
Mott insulator state, the ground state must break trans- 
lation symmetry to enlarge the unit cell such that the 
number of electrons per unit cell becomes even. Here we 
present a counter-example. The ground state of Sp{2N) 
when is large can be featureless Mott insulator with 
odd numbers of electrons per unit cell. 

We want to introduce a large-A/' limit and maintain 
the SU (2) gauge structure. The simplest way to do this 
is to introduce A^ flavors of fermions (later in this paper 
we also denote A^ as Nf to emphasize it represents the 
number of flavors of fermions): 



JxX 

hi 



l,2,---N. 



Then the Lagrangian of the A^-flavor model is: 
L = ^Vr(iSt-4,,r')V? 



13 4'^»J 

(ij) 



'-TvUijUj. + {4''UijiP] + h.c.) 



where the repeated index a is summed. The large-A?' 

model is clearly invariant under the SU (2) gauge trans- 
formation. Comparing with the original model Eq. (7), 
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after integrating out the fermions, the effect of the TV- 
flavor is to put a factor of N in front of the boson fields 
Lagrangian, which controls the fluctuation of them when 
A'^ is large. 

Now we want to find out the corresponding spin model 
for this large- -/V theory. It should be some generalization 
of the spin-1/2 Heisenberg model. To do so, we need to 
integrate out the boson fields. If we go back to the / 
picture, we have: 

i 



{ij) 



where 



X 



Here one immediately see the constraints become Xii = 
N and r??? = 0, or 



j-of fa 



N 



Jia ct0Ji 



at 
03 



h.c. = 



(15) 



for each site. The Hamiltonian under these constraints 
are: 



(ij) 



E 

(ij) 



'an 



{ _ f of fb fdf fb 
y JiaJil3Jjl3Jja 



-itliy^fu + f:lfu;if^p + n) h?) 

We want to understand the symmetry of this model 
and try to rewrite it in a form of spin coupling. Here, 
the symmetry that we are studying is not the local gauge 
invariance, but some global physical symmetry in analogy 
with the spin rotation symmetry for the A'^ = 1 model. 
Motivated by Eq. (11-13), we construct all the gauge in- 
variant bilinear forms of for each site (neglecting the 
site label). Let 



ta 



ft 
-f\ 



at 



The SU{2) gauge invariant bilinears are 
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What is the group generated by these S operators? First 
let us count how many of them there are. For 5"''+ or 
S"***", the label is symmetric for a and 6, so there are 
^^^^"'"■^ operators of each type. For S""^^ , the labels are 
not symmetric, so there are simply iV'^ of them. Totally 
we have iV(JV + 1) + A''^ = 2N'^ + N oi them. One can 
further examine their commutation relations: 



1^ gabs ^cd3 j 



^^bc i^adS ^adgcbS'^ 

ac i^bdS _|_ ^ad gbcS 
_|_^6c i^adS _|_ ^bd igacS^ 



0, 

1 

2 
1 

2 

_ 1 (^gadgbc 
1 



) 



2 ^ 



These are the relations for SP{2N) algebra, so all the 
S operators arc the 27V^ + N generators which generate 
an SP{2N) group. When N = 1, SP{2) is isomorphic 
to 50(3). After some rearrangements, the Hamiltonian 
Eq. (17) can be rewritten as: 



H 



N 



{ij) 



^ Qab-\- Qab- 

2 i ^ 



Qab— Qab-\- 

2 i ^ 



ab3 Qba3 



(18) 



^ = - E il {^t;'v^ + n;'x^ (16) If we define: 



^ab- 



gabl _ gbal _ _ ^^^^+ _|_ 

Qba2 ^ / Qab-^ Qab—\ 



gab2 

then the vector 

Qdb ( oabl Qab2 cab3\ 

'-'i ~ I'-'i ' '^i ) 

can simplify the Hamiltonian to: 



H. 



Sp{2N) 



E ^ij cab Q 



ba 
3 



(19) 



Here we have to mention that the three components of 
Sf" are actually not on the same footing, since the first 
two are symmetric with respect to the flavor labels but 
the third one is not. It is a simple task to check that 
all the generators commute with the Hamiltonian. And 
one can even check that the SP{2N) is indeed the full 
symmetry group of it. 

Now we want to know what the physical Hilbert space 
is. It should be the subspace of the enlarged Hilbert space 
in which the constraints Eq. (15) arc satisfied. When 
A^ = 1, the physical Hilbert space is the spin-1/2 Heisen- 
berg's, where we have only two states on each site: 

IT), li) 
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When N = 2,we have 5 states on each site: 
|1t2t), |lT2i), \li2^), 
-^(|ln2o)-|lo2n)) 

Here, for example, |1|2|) represents a ferniion of flavor 
1 and spin up and another fermion of flavor 2 and spin 
down; and |l||2o) represents a fermion of flavor 1 and 
spin up, another fermion of flavor 1 and spin down, and 
no fermion of flavor 2. One can check the dimension of 
the Hilbert space on each site for larger N's: 

N=3 dimension = 14, 

N=4 dimension = 43, 

N=5 dimension = 142, 

N=6 dimension = 429, • • • 

This Hilbert space turns out to be an irreducible repre- 
sentation of the SP{2N) symmetry group. If wc label an 
irreducible representation by its highest weight state for 
a particular Cartan basis. The Cartan basis for SP{2N) 
can be chosen to be the z-component spins of each flavor: 

S""^, where a = 1,2,- •• ,N, 

with N generators. Then the highest weight state in our 
Hilbert space is simply: 

|l-|-2-|-3-|- . . . iV|) 

IV. PHASE DIAGRAM OF THE Sp{2N) MODEL 
ON 2D SQUARE LATTICE 

In this section we would like to calculate the zero tem- 
perature phase diagram for the Sp{2N) model Eq. (19) 

on 2D square lattice. Wc assume that only the nearest- 
neighbor couplings and the next-nearest-neighbor cou- 
plings are non-zero 

Ji,i+x — Ji,i+y — Jli 
Ji,i+x+y — Ji,i+x—y — ^2- 

We also assume Ji + J2 — 1- 

There are two mean-fleld approaches to the Sp{2N) 
model. In the first mean-field approach, we use the 
ground state )) of a trial Hamiltonian 

i?tHa. = E^mf -Sf (20) 

as the trivial wave function and obtain the mean-field 
ground state of the Sp{2N) model Eq. (19) by minimizing 

as we vary the variational parameters m"^. The ob- 
tained ground state corresponds to a Sp{2N) spin polar- 
ized state. One can show that the ground state energy 



obtained in this mean-field approach is always of order 1 
per site in the large N limit. 

In the second mean-field approach (the projective 
construction), [1, 25] we start with a fermion trial Hamil- 
tonian 

Htrial = ^^fuij^P] (21) 
(ij) 

where Uij = u^^t^ are two by two complex matrices that 
satisfy 

Aj = '"ii' "ii = imaginary, u\j = |i=i,2,3real 

Let |$TOeara) be the ground state of the above trial Hamil- 
tonian. u\^ = are chosen such that the constraints 

ipi^T^tpi = are satisfied on average: 

U/.of i,/b|,f,("«) \ _ Q 

We then project \^mean) to the physical Hilbert space 
and obtain P|$meara)- The physical Hilbert space is 
formed by states \phys) that satisfy the constraints 
Eq. (15) or 

The projected wave function P\^mean) is our trial wave 
function with Uij as variational parameters. Minimzing 

{^tei'„\PHsp^2N)P\^teii) 

by varying mj , we obtain the approximated ground state. 

Since in the large- TV limit, the fluctuations of the La- 
grangian multiplier Qq ^{1) are weak, we expect that re- 
moving the projection P only causes an error of order 
1/N. Also, other mean-fleld fluctuations are weak in the 
large-TV limit, so we expect that the minimized mean- 
field energy 

E = {^meaL\Hsp(2N)\^meim) 




to be the true ground state energy at the leading order 
in the large- A?^ expansion. Here 

Xij — \^mean\Xij\^^eo.n/i 'lij — \^ mean \ ' lij \^ mean/ 

note that the above minimized energy is of order —N 
per site. Thus the states obtained in the first mean-field 
approximation Eq. 20 cannot be the ground state. 

Within the SU{2) projective construction and using 
the translation invariant ansatz Uij with only ncarest- 
and next-nearest-neighbor couplings, we find many local 
minima of -EnjMfvt^j + xtfxfj) as we vary mj. 
We plot those minima in Fig. 1 as functions of J2 (note 
Ji + J2 = 1). 
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As we change J2, the energy of the state changes 
smoothly along each curve. So there is no quantum phase 
transition as we move from one point of a curve to an- 
other point on the same curve. The ansatzs on the same 
curve belong to the same phase. However, if two curves 
cross each other, the crossing point represents a quantum 
phase transition. This is because the ground state energy 
is not analytic at the crossing point. If the slopes of the 
curves at the crossing point are different, the quantum 
phase transition is first order. If the slopes at the cross- 
ing point are the same, the quantum phase transition is 
second order. 

From Fig. 1, we see second-order (or continuous) phase 
transitions (at mean-field level) between the following 
pairs of phases: (A,D), (A,G), (B,G), (C,E), and (B,H). 
We used to believe that all the second-order phase tran- 
sitions are caused by symmetry breaking. So a natural 
question is what symmetries are broken for the above five 
second-order phase transitions? 

It turns out that, except phase (D) and phase (E), all 
other phases have the same symmetry. In other words, 
the projected ground state wave functions P\^mean) for 
the ansatz Uij associated with those phases have identical 
symmetry. Thus the three continuous transitions (B,G), 
(B,H) and (A,G) do not change any symmetries. It was 
pointed out that those phases, despite having the same 
symmetry, contain different quantum orders [12]. The 
projective symmetry group (PSG), defined as the invari- 
ant group of the ansatz Uij , is introduced to describe this 
new class of orders [12]. 

The ansatzs on the same curve have the same PSG 
and correspond to the same quantum phase. On the 
other hand, the ansatzs on the different curves have dif- 
ferent PSG's. We see that a quantum phase transition is 
characterized by a change in PSG. Those quantum phase 
transitions represent a new class of phase transitions be- 
yond the Landau's symmetry breaking theory. Other ex- 
amples of phase transitions beyond the Landau's theory 
can be found in Ref. [19-23, 26]. 

In the following we will discuss quantum orders (or 
PSG's) for mean-field phases in Fig. 1. The PSG's for 
those quantum orders are labeled by labels which look 
like Z2A0013 [12]. 

The phase (A) [27] is the vr-flux state, or the SU2BnO 
state 

=iXi Ui^i+y =i(-)*"x 4 0- (22) 

The low energy excitations are described by massless 
Dirac fermions with a linear dispersion and gapless SU (2) 
gauge fluctuations. Therefore we also call such a state 
S'C/(2)-linear state. 

The phase (B) [12] is a state with two independent uni- 
form RVB states [1] on the diagonal links. The gapless 
fermions have finite Fermi surfaces. The fermions inter- 
act with SU{2) X SU{2) gauge fluctuations. Such a state 
is called SU{2) x 5t/(2)-gapless state {SU{2) x SU{2) 
indicates the low energy gauge group and "gapless" indi- 




FIG. 1; The mean-field energies for various phases in a Ji- 
J2 spin system. (A) the yr-flux state (the S'C/(2)-linear state 
SU2BnO). (B) the SU{2) x SU (2) -gapless state. (C) the 
SU{2) X SU {2)-\mear state. (D) the chiral spin state (an 
5'(/(2)-gapped state). (E) the [/(l)-Unear state Eq. (26) which 
breaks 90° rotation symmetry. (F) the L''(l)-gapped state 
UlCnOOi. (G) the Z2-linear state Z2Az«13. (H) the Z2-linear 
state Z2A0013. (I) the uniform RVB state (the S'(7(2)-gapless 
state SU2AnO). 

cates finite Fermi surface). Its ansatz is given by 

Ui.i+x+y =XT^, Ui^i+a:-y =XT^ , «[) =0- (23) 

The phase (C) [12] is a state with two independent tt- 
flux states on the diagonal links. It has SU{2) x SU{2) 
gauge fluctuations at low energies and will be called an 
SU{2) X S'[/(2)-linear state. Its ansatz is given by 

Ui,i+a:+y =x{t^ + T^) Ui^i+^-y =x(t^ - T^) 

4 =0 (24) 

The low energy excitations are SU{2) x SU (2) gauge fluc- 
tuations and massless Dirac fermions. 
The phase (D) is the chiral spin state[14] 

XiA+x=iXi, XiA+y ^■iXii-T"' , 4=0' 

XiA+j:+y = - «X2(-)*", XiA+a:-y =iX2(-)'". (25) 

Both fermionic excitations and SU{2) gauge excitations 
are gapped. The gap of the SU{2) gauge excitations is 
due to an SU (2) Chern-Simons term. 

The phase (E) [12] is described by an ansatz 

Ui,i+x+y =XlT^ + X2T^ Ui^i+x-y ^XlT^ " X2T^ 

UiA+y ^rjT^ flo =0 (26) 

which breaks the 90° rotation symmetry. It is a U{1)- 
linear state, i.e., the low lying excitations are massless 
C/(l) gauge fluctuations interacting with massless Dirac 
fermions. 
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The phase (F) [12] is described by the following 
UlCnOOx ansatz 



3 

al =A, 



'.,i+x—y — X''" 



,1,2 







(27) 



The UlCnOOa; state can be a C/(l)-linear state where 
fermions are gapless with a linear dispersion relation (if 
Oq is small) or a [/(l)-gapped state where the fermions 
are gapped (if Og is large). The state for phase (F) turns 
out to be a [/(l)-gapped state. The only low energy ex- 
citations are massless U{1) gauge bosons. 

The phase (G) [12] is described by the Z2Azzl3 ansatz 



UiS+y ^XT^ +VT^, 



^i^i—x+y t 
1,2.3 



=0. 



(28) 



The SU{2) gauge structure is broken down to a Z2 
gauge structure. Hence there is no gapless gauge fluc- 
tuations. The only low energy excitations are massless 
Dirac fermions. Such a state is called Z2-linear state. 
The phase (H) [18] is described by the Z2A0013 ansatz 



4 ^ 0, aY = 0, 

UiA+x = XT^ + VT^, 
UiA+x+y = +1T^, 



Ui.i+y = XT^ - VT^, 
Li-a^+y = +1T^- (29) 



It is also a Z2-linear state. 

The phase (I) is the uniform RVB state [1], i.e., the 
S';7(2)-gapless state SU2AnO 



aj, = 0. 



(30) 



It has gapless SU{2) gauge fluctuations and gapless 
fermionic excitations that form a finite Fermi surface. 

From Fig. 1, we see continuous phase transitions (at 
mean-field level) between the following pairs of phases: 
(A,D), (A,G), (B,G), (C,E), and (B,H). For the three 
continuous transitions (B,G), (B,II) and (A,G) that do 
not change any symmetries, we observe that the SU (2) 
gauge structure in the phase (A) breaks down to Z2 in 
the continuous transition from the phase (A) to the phase 
(G). The SU{2) x 5/7(2) gauge structure in the phase 
(B) breaks down to Z2 in the two transitions (B,G) and 
(B,H). 



V. A STABLE ALGEBRAIC SPIN LIQUID - 
S(/(2)-LINEAR SPIN LIQUID 

In this section, we will study the iS'[/(2)-linear state 
given by Eq.(22) which describes the phase A in phase 
diagram Fig. 1. In the mean-field theory, the interactions 
between the excitations are ignored. In this section, we 
include those interactions and study how those interac- 
tions affect the low energy properties of the mean-field 



X 



1 



FIG. 2: We choose the unit cell of S'(7(2)-linear state to con- 
tain jouT lattice sites when we go to the fc-space. 



S'?7(2)-linear state. We will show that, after including 
those interactions, the gapless excitations in the mean- 
field 5J7(2)-linear state remain gapless, which leads to a 
stable algebraic spin liquid. 



A. The low energy effective theory of the 
Sf/(2)-linear state 

To obtain the low energy effective theory of the SU (2)- 
linear state in the continuum limit, we choose the unit 
cell as in Fig. 2. It contains 4 sites; each site has spin up 
and down, so totally 8 fermions. Let us write down the 
mean field Hamiltonian: 



Pi-\-x,2 



h.c. 



k 

' iX [V'LV'fe2 - 4iV'fe3 + ipU'ipkoe'^" - V'LV'fcie''''' 



h.c. 




(31) 



where 



/ 


i 


— ie '■'^^ 


i — ie 









+ ie^^- 





- i 


+ ie' 


iky 


—i 


+ ie'^^y 





i 


— ie" 




V 


i 


— ie'*^" 


— i + ie'-'^^ 








(32) 



Here we have assumed the lattice constant to be 1/2, so 
we have — tt < kx,ky < tt. Note that here i/'fc.i is actually 
S't/(2)-doublet, corresponding to the spin up and down 
components in the /-formalism; so still totally 8 fermions. 
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After some rearrangements: 
M = 2x- 



/ 

— e 

—e^'^y sinfe 



sin kx — e » sin fc 



sin fca; 







V 







e^''y sin k,. 






nikx 



sinfcx 





sin /cj, 




(33) 



In the continuous limit, the energy spectrum for fermion 
is characterized by a single fermi point at (0,0). When 
w 0: 



energy fluctuations. Understanding those terms is the 
key to understand the low energy behavior of the model. 
Those terms must be consistent with the underlying lat- 
tice symmetry. So in the following, we will study the 
symmetry properties of the effective theory Eq. (40). 
We will show that none of terms allowed by the sym- 
metry are relevant at low energies. None of those terms 
can cause infrared instability. As a result the mean-field 
SU (2)-linear state leads to a stable algebraic spin liquid. 



B. Space translation and rotation symmetry 



M = 2x 




(34) 



We can do an extra rotation to make it the usual form 
of Dirac fermion: 



R 



71 tI 
71 71 

V 71 71 y 



(35) 
(36) 



Then 



Hmean = '0'^7O [ikxjl + ikyj2] 1p 



(37) 



where the 7 matrices in Euclidean space are: 



70 



■ . 



-0-3 



71 



-CTl 



72 



-<J2 

(38) 



Here we can also introduce the other two 7 matrices: 



73 = 



75 = 70717273 = i 



-I 



(39) 



where / is the 2 by 2 identity matrix. Notice that both 73 
and 75 anticommute with all space-time components of 
7 matrices: 70, 71 and 72. We should also include gauge 
field fluctuations above the mean-field theory. The full 
lagrangian is: 



1 



L = i;{d,- zay) + ^Tr [/^^/^J 



(40) 



where ijj = ip'^'fo- So the low energy effective theory of 
the S'i7(2)-linear state is a QCD3 with a SU(2) gauge 
field and two massless 4-component Dirac fermions (or 
^Nf 4-componcnt Dirac fermions which form Nj SU{2) 
gauge doublets in the large Nf limit, notice iV/ = 1 is 
the physical Heisenberg model's case). 

The • • • in Eq. (40) represents other terms which may 
be generated by the interactions as we integrate out high 



Now let us think about the corresponding lattice sym- 
metry in continuous limit. Firstly let us discuss transla- 
tion by one lattice site along x-direction T^, in terms of 
the lattice fields: 



:{ 



(41) 



It seems that the translation symmetry is broken since 



Uij is not invariant: 



{ 



= -ii-Y'^x 



(42) 



But as shown in Fig. 3 one can do an extra local SU{2) 
gauge transformation Wt^ to transform Uij back. Let 



Then 



W- 



t 



where 



II 

i,i+y 



= W 



-u' 



= i{-y''X = Ui,i+y 



(43) 



(44) 



(45) 



Here wc point out that the combination of Wt,, and T,,: 
is a transformation leaving Uij invariant. We call such a 
transformation an clement of PSG. PSG, by definition, 
is the collection of all transformations leaving the ansatz 



invariant. Here in S'C/(2)-linear state, ^ is zero, 
so it is also invariant. 

After choosing the unit cell as in Fig. 2, in terms of the 
four component ip fermion, the combination of Wt^ and 
Tx transforms: 







/O 1 







_ 


1 





V'^2 














\0 


-1 



Ipil 

V'i2 



(46) 



Here wc assumed in the continuous limit, ipi and V-'i+x 
are on the same position. We can do the extra rotation 
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A 

























(a) 



(b) 











1 


1 


1 




-1 


-1 


-1 




1 


1 


1 





(c) 



(d) 



FIG. 3: Figures (a)-(d) illustrate how the SC/(2)-linear state 
respects translation symmetry Tx'. (a) original ansatz (b) after 
Tx (c)do Wt^ (d) go back to original ansatz. 



to transform into ip, the usual Dirac fermion: 




(47) 



Eventually we know that the PSG element Wt^ o T^, 
which is the lattice symmetry, if translated into contin- 
uous limit, is the internal symmetry —73. Note that the 
minus sign actually corresponds to a global SU(2) trans- 
formation: Wi = — 1. So we have the correspondence: 



73 



(48) 



Similarly one can study the translation by one lattice site 
along y-direction; we find: 



Wt^ o Ty < — > 75 
For the reflection : x ^ —x: 

Wp^ o < — > 71 
For the reflection Py : y ^ —y: 

Wp^ o Py < — > 72 



(49) 



(50) 



(51) 



For the reflection P^y : x 



xy 



y,y -> x: 



2(71 -72) (73 +75) (52) 



Note that Tx,Ty, P^, Py, Pxy already give the full space- 
time symmetry. For example, rotation by 90 degree i?go : 



x^y,y 



is a combination of Pr and Pn 



xy 



R90 — Pxy ° Px 



C. Time-reversal symmetry 



(53) 



Now let us study the time-reversal symmetry T . In 
terms of spin operator: 



T : S 



(54) 



One should be cautious that T is not a usual linear op- 
erator, instead it is an anti-linear operator: 

Ti = -iT (55) 

What does this operator correspond to in terms of lattice 
fermion and bond variable Uijl We know that 



1 



(56) 



So the corresponding transformation on / fermion is: 

T:f^ia2f, T : /t ^ /t(^a2)^ (57) 

= P{-icr2)^ ('7l,-CT2,cr3) («cr2)/ 
= P (-CTl,-cr2,-Cr3)/ 

= -S (58) 

Here we used the anti-linear property of T operator. 
Therefore the T transformation on tp fermion is: 



T : V 



fl 

-n 



(59) 



Our convention of notation is that in terms of / fermion, 
we use a to denote Pauli matrices; while in terms of ^ 
fermion, we use r to denote them. 

What is the time-reversal transformation on Uij? Here 
we notice that Uij has two meanings: in Eq.(7) it means 
the operator defined as in Eq.(6); while in Eq.(14) it 
means the average value of the operator on mean-field 
ground state. Let us take a look at how the operator 
transform under time-reversal: 

^fK^<J2y{^<J2)fj = fUj=Xij (60) 

(61) 
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Therefore we know that the JJij, as an operator, is in- 
variant under T. Here it is helpful to write down tjij in 
terms of V' operators: 



Xij Vij 

4j -x^ 



-V^iVl - {ir2rj){ir2ri)^ (62) 



Notice Eq.(59), under T transformation, the first term 
transforms into the second term, and the second trans- 
forms into the first. So the whole Uij is invariant. One 
can also check that together with transformation Eq.(59), 
the Lagrangian Eq.(7) is invariant under T. This is ex- 
pected since the original Heisenberg Hamiltonian Eq.(l) 
is T invariant. 

But things arc different if Uij condense, more specif- 
ically, if it has some non-zero average value: Uij = 
{'^lUijl^). Because T is an anti-linear operator, it trans- 
forms Uij into U^j . Wc know that for two arbitrary states 
\a), 1/3), the anti-linear property of T gives 



{a\/3) = {Tp\Ta) 
Therefore for an arbitrary linear operator O 

{a\d\(3) = {6^a\P) = (Tf3\Td^a 



(63) 



{Tp\TO'iT-'^\Ta) 
(64) 



let d = Xij: and \a) = |/3) = 1^*), one immediately see 
that under T transformation, the average value of Xij 
operator transforms as: 



Xij = (*IXiil*) 



X'ij = {T^\Xij\T^) 

= mlm = xtj 



(65) 



Here we know that the mean- field variable Uij, as the 
average value of Uij operator, transforms into U^j under 



T. This is consistent with our understanding of anti- 
linear T operator, since it transforms any c- number into 
its complex conjugate. 



T : Uij ^ TUijT-' = Utj 



In general, Uij can be written as: 



Uij 



u^ii-ro 



(66) 



(67) 



where tq is identity matrix. For spin rotation invariant 
system, one can show that u^j is pure imaginary, while 
uij i = 1, 2, 3 are pure real. 

How does mean-field Hamiltonian transform under T? 
Here we need combine transformations on tjj and Uij : 

V'l Uij Ipj ^Ipf (iT2 ) Uij {iT2)lpj 

(using Eq.(67)) =i;f{-Uij)rj 

=i^l{Uijfi>i 



(Uij = uU 



'-i^^jU^jii^i 



(68) 



Comparing with the term in original mean-field Hamilto- 
nian ipjUjiipi, one concludes that the transformation of 
mean-field hamiltonian under T can be simply expressed 



as Ui- 



U't^ , with no fermion transformation. If Uf^ and 



original Uij can be related by a SU (2) gauge transforma- 
tion, the system has time-reversal symmetry; otherwise 
the T was broken. 

For our S'?7(2)-linear state, U^- and Uij are indeed re- 
lated by a gauge transformation Wt- 



Wt : Wi 



-(-) 



(69) 



After choosing unit cell, in terms of four-component 
fermion ■)/;, 



Wt 



In terms of usual Dirac fermion ip: 
1 











II 








1 V'i2 









(70) 



-1 



-1 



V' = 7oV' 



(71) 



1 



The combined transformation Wt ° T leaves Uij invari- 
ant. So Wt o T is another element if PSG. And we know 
in terms of Dirac fermion ip, 



Wt oT : ip ^ tj/ = (2x2)70?/' 



(72) 



This is very similar to the time-reversal transformation 
in usual Dirac field theory. 



D. Spin rotation: "charge conjugation" 

There is another important symmetry, the global spin 
rotation. Let us think about rotation around y-axis by 
180 degree: 



In terms of / spinon. 



Rs 



/t 



(73) 



(74) 



In terms of ip fermion. 



R. 



■spin 



fl 
-fl 



t = {mW (75) 



It seems Eq.(75) is identical to Eq.(59), which means T 
and Rspin are identical. This is obviously wrong, the 
difference here lies in the fact that T is anti-linear but 
Rspin is linear. 
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In terms of 4-component Dirac fermion ijj: 

Rspin : = R^V ^R^(ir2)V* = (ir2)RtR>* 

={iT2){-7i75)r (76) 

The transformation rule above is in real space. In mo- 
mentum space, 



Rs 



-ikj 



^ (-ir2)(7i75)-^X]^"''''^J* 
j 



(77) 



i.e., the spin rotation transformation -Rspin also flip the 
sign of momentum! 

Now think about how Uij transforms under Rspin ■ The 
form of Eq.(62) is explicitly Rspin invariant, so Uij is 
R 



spin 



invariant. This is expected since wo have spin ro- 
tation invariant x and f}. This is actually the point of this 
slave-boson mean-field approach. We choose the mean- 
field variables to be spin rotation singlets, such that even 
it acquires non-zero average value, it doesn't break spin- 
rotation symmetry, which describes spin liquid states. 
Moreover, since Rspin is a linear operator, we know that 
the average value of Uij is also Rspin invariant. 



Rspin • Uij 



Rspin Uij Rgpin 



(78) 



This is consistent with the usual property of linear oper- 
ator: it commutes with c-numbers. 

What is the transformation on mean-field Hamilto- 
nian? We should put transformations on fermion Eq.(75) 
and Uij Eq.(78) together: 

Rspin ■ i>\Uiji^j ^IpJ {iT2)^ Uij{iT2)ljjj 



using Eq.(67) 



--i'l{-u:j)rj 

--i^]uyi = ^\Uji^i 



So Rspin is also an element of 



This is just another term in the original mean- field 
Hamiltonian. We conclude that Hmean is invariant under 

Rspin ■ This again is expected since onr spin- liquid states 
should not break Rg 
PSG. 

Here we point OTit that the form of Rspin transfor- 
mation Eq.(76) is very similar to the charge-conjugation 
transformation C in usual Dirac field theory. In fact we 
can denote Rspin as C in our later discussion. We sum- 
marize the transformations on Dirac fermion in Table I. 



PSG element 


Transformation on fermion 


^'^ = Wt^ 


73 




78 




71 




72 




1(71 -72)(73 + 75) 


T^^G = WtoT 


— » {iT2)fo'4>* (anti-linear) 




tp — > —{iT2)fi^5ip* (linear) 



TABLE I: Transformations of Dirac fermion ip, here C is Rspin 
and Wc is identity. 



be generated after including quantum fluctuations? If 
yes then the mean-field theory cannot describe the real 
physical system at all since the low energy behavior is 
completely different. In this section we will show that 
this is not the case. 

We know that for a certain mean-field state, it has a 
certain PSG, which is the symmetry of the mean-field 
hamiltonian. After we include fiuctuation, this symme- 
try will still be respected. Suppose we go through an 
renormalization group process to find the low energy the- 
ory, any counter-term explicitly breaking PSG will not be 
generated when taking care of fiuctuations. 

This can be viewed in the following way: let 
L{ipi,Uij,aQj) in Eq.(7) be the Lagrangian describing 
the dynamics of fermion and gauge fields. In the orig- 
inal theory, we have a huge "symmetry" group leaving 
L invariant, for example, the translation along ar-axis by 
one lattice site Tr: 



L{^i, Uij) = L{T,yjiT-\T,UijT-^) 



(79) 



or an arbitrary local SU (2) gauge "symmetry" transfor- 
mation W (here the meaning of quotion mark is that 

gauge "symmetry" is not a physical symmetry, instead it 
is just a many-to-one bad labelling.): 



L{^pi, Uij) = L{WixPi, WiUijWh 



(80) 



But after Uij condense, things are different. The above 
huge "symmetry" group will be "spontaneously" broken 
(the unbroken state must have Uij = 0). PSG is the re- 
maining unbroken "symmetry" group after this symme- 
try breaking. PSG is defined as the collection of all trans- 
formations leaving the Lagrangian L{'ipi, Uij, Qq^) invari- 
ant and also leaving average value Uij invariant. Let P 
be an element of PSG, then 



E. Transformations of fermion bilinears 

Let us focus on SU (2)-linear state. We know that the 

mean-field Hamiltonian is characterized by a SU (2) dou- 
blet of massless Dirac fermions. But can the mass gap 



LitpijUij) 
Uii 



L{Pi,iP-\PU 
PUijP-^ 



1-3 



(81) 
(82) 



We can consider fiuctuation around the average value 
Uii- 



u. 



1-3 



u 



1-3 



SU, 



1-3 



(83) 
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Plugging into Eq.(81) 

L{iPi, Uij + 6Uij) = L{Pi)iP-\ P{Uij + SUij)P-^) 

= L{P^iP-\Uij + PSUijP-^) (84) 

So the Lagrangian for the fluctuations Lijpi, 6Uij) must 
be invariant under PSG transformations. 

We want to find out the transformations of fermion 
bilinears under PSG. Because if they all transform non- 
trivially under PSG, they will be all forbidden. That is 
why the fermions remain to be massless after including 
fluctuations. 

Here wc will consider only the fermion bilinears of form 
The bilinears of forms tpip and ijj^tfj^ are not invari- 
ant under the spin Sz rotation and are not allowed in the 
effective Lagrangian. Since we are using 4-component 
fermion, there should be 4 x 4 = 16 different fermion bi- 
linears of form tp^tp. From Table I, it is quite easy to flnd 
the transformations of all fermion bilinears under PSG. It 
turns out that among these 16 bilinears, as shown in Ta- 
ble II, there are four 1-dimensional representations and 
six 2-dimensional representations of PSG. All the fermion 
bilinears transform non-trivially under the PSG. So per- 
turbatively, the fermions remain massless after inclusion 
of fluctuations. 

Now let us consider the fermion bi-linear terms that 
also contain a single spatial derivative. Those terms rep- 
resent marginal perturbations when Nf = oo. From the 
table II, we see that the only term that is allowed by 
the PSG is t/jdx'fiip + ipdy'-f2'>P- AH other terms are for- 
biden. The reason is as follows: The 1-dimensional rep- 
resentations together with a spatial derivative cannot be 
Lorentz singlet, so are ruled out. Among 2-dimensional 
representati5)ns together with a spatial derivative, only 
^dx^itp + 4>dy'y24' is invariant under translation T^^'^ 
(in fact, it is invariant under full PSG); all others are 

ruled out. But the term ipdi'yitp + tpd2j2''P, which is al- 
ready present in Eq.(40), only changes the velocity of 
the fermions. In a RG study, this counter-term means a 
wavefunction renormalization. The low-energy effective 
theory Eq.(40) remains valid. 

The next question is, will there be 4-fermion interac- 
tion terms? The answer is yes. For example, we choose 
a certain 1-dimension representation in Table II, say tpip, 
then couple this term to itself to make a 4-fermion in- 
teraction. It is obvious that this 4-fermion term is PSG 
invariant, wliieli is allowed in the Lagrangian. Will this 
kind of 4-fermion term change the low energy behavior 
drastically? The answer is no. This is because we are 
in 2-1-1 space-time dimension, and by power counting 4- 
fermion terms are of dimension 4, so they are irrelevant 
couplings. Same argument can be done for fermion bilin- 
ear terms with second order derivatives; those terms are 
also irrelevant. 

In summary, we have discussed the possible fermion 
self- interactions. Our conclusion is that in perturbative 
sense, these fermion self-interactions will not change the 



low energy behavior from the mean- field result. We may 
say that the SU (2)-linear mean-field state is stable under 
fermion self-interactions. 



F. Emergent 5p(4) physical symmetry 

In this section we will discuss the emergent symmetry 
for SU (2)-linear phase whose low energy effective theory 
is Eq.(40). We already know that there are two fermion- 
chiral symmetry generators 73 and 75, and they are anti- 
commuting. (Here please note that we are talking about 
fermion-chiral symmetry, which is different from the spin- 
chirality in the later discussion about chiral spin liquid.) 
Therefore the symmetry of SU (2)-linear phase contains 
at least a global SU{2) fermion-chiral Lie group whose 
generators are 73,75 and ^7375. 

We also know that the theory should be global SU (2) 
spin-rotation invariant, since we are talking about a spin 
liquid phase here. Thus there should be at least another 
SU{2) spin-rotation symmetry group. However after the 
particle- hole transformation we made in Eq.(8), this spin- 
rotation symmetry was hidden in our formalism. 

The full physical symmetry group of 5?7(2)-linear 
phase should contain both the SU (2) fermion-chiral and 
SU{2) spin-rotation as its subgroups. The naive guess 
for the fufl group is SU{2) x SU{2) but it turns out to 
be wrong. We will show in this section that the cor- 
rect full symmetry group is 5^(4). We noticed that the 
same 5*^(4) symmetry was found earlier by Tanaka and 
Hu[28] by viewing the 7r-flux state as a fermionic mean 
field state, i.e., ignoring the effect of SU{2) gauge field. 
Here we clarified the gauge field effect and obtained the 
same global flavor symmetry. Then we can classify all 
the fermion bilinears according to their transformation 
rules under 5^(4) group. 

Later, we will show that after including the SU{2) 
gauge fluctuations, the 5'?7(2)-linear state remains gap- 
less and the correlations between various operators re- 
main algebraic. But the exponents of algebraic correla- 
tions may be modified by the SU{2) gauge interaction. 
Classifying fermion bilinears according to their transfor- 
mation under the 5*^(4) group is very important in under- 
standing the scaling properties of those operators. The 
operators that belong to the same irreducible Sp{4) rep- 
resentation will have the same scaling dimension. 

First we consider the spin rotation group; basically it 
will mix tp and ■0* . To make the spin rotation trans- 
formation explicit, it is convenient to reintroduce the tp 
formalism in Eq.(lO): 

i, = ia2r = /jt ) (85) 

Let us look at lattice fermion at certain site. If we put 
ip (2-component, corresponding to spin up and down) and 



13 



nPSG 



+ 1 



+ 1 



+ 1 



-1 



+ 1 



+ 1 



+ 1 



yPSG 



-1 



+ 1 



-1 



+ 1 



+ 1 



jPSG 



-1 



+1 



-1 



+ 1 



+ 1 



yPSG 



+ 1 



+ 1 



-1 



+ 1 



+ 1 



T 



PSG 



C 



PSG 



-1 



-1 



V'7072'!/' 



-1 



+1 



-1 



+1 



-1 



+1 



+1 



+1 



+1 



+1 



+1 



+ 1 



+1 



-1 



-1 



-1 



+ 1 



#7173^' 

#7275^/' 



+ 1 



+ 1 



+1 



+ 1 



i'07i75i/) 

#7273^' 



+ 1 



+1 



+1 



-1 



-1 



+1 



-1 



TABLE II: Transformations of 16 fermion bilinears. All transform non-trivially. 



(2-component) together to form a 4-component vector: 

then it is straightforward to write down the spin rotation 
transformation. For example, the rotation around z-axis: 



thus 



V'- 



Therefore 



9/2 



ie/2 Q 







-ie/2 



(88) 



where the identity matrix labels the internal space of 
(spin up and down), while e^^'^^l'^ acts on the space 

mixing i/; and "p. 

What about rotation along t/-axis? Suppose we do an 

infinitesimal transformation 



it implies: 



h 

h 

V'l ■ 

V'2 ■ 



/t + 



V'2 - ^v-r 



(89) 



(90) 



in terms of ^': 



(87) 

Similarly the rotation along a;-axis is: 



(91) 



(92) 



(93) 



To summarize, we know that the spin rotation is acting 
on the space mixing i/) and i]). 

Let us go to continuous limit, and consider the 4- 
componcnt Dirac fermion SU(2) doublet in Eq.(36). 

Again we write it together with t/) = (—7175)^0-2^/'*, 
where the (—7175) is inserted to make the spin-rotation 
have a simple form: 




(94) 



Note that actually ^ has 16 components and 16 = 4x2x2 

where 4 is the number of Dirac components, the first 2 
is for SU{2) gauge doublet and the second 2 is for the 

space mixing ^ and ^. From the above discussion, the 
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space mixing ip and ip is actually spin rotation space. A 
generic transformation G on fermion field can be written 
as a transformation on ^: 

irac ^ Ggauge G spin {^^) 

where the transformations with subscripts act on each 
corresponding space. 

The three spin rotation generators are, from above dis- 
cussion: 

10 10C71, 1(g) 10 (72, 10 10(73 (96) 

The fermion- chiral generator 73 is acting on ijj. One 

can easily check that while acting on since -ip = 
(— 7i75)«cr2^*, the generator has the form: 73 1 
£73. Similarly one can find the other two generators of 
fermion- chiral transformation. In summary, they are: 

73 10(73, 75 100-3, ^7375®!®! (97) 

Now if we do commutations between Eq.(96) and 
Eq.(97), the full set of symmetry generators can be found: 

1 1 cr, 73 1 cr, 75 1 cr, ^7375 10 1 

(98) 

Totally 3 + 3 + 3 + 1 = 10 elements, which satisfy 5^(4) 
algebraic relation. 

Here one thing we need to mention is that the three 
gauge transformation generators: 

1 T 1 (99) 

will also keep the Lagrangian Eq.(40) invariant. But they 
arc gauge transformations and should not be taken as 
physical symmetries. 

We just showed the emergent Sp{4:) global symme- 
try. Can the emergent continuous symmetry group larger 
than 5*^(4)? The answer is no, as one can see in the fol- 
lowing. Wc have totally 8 components of fermions, and 
they form four SU{2) gauge doublets. For global sym- 
metry we should only consider transformations invariant 
in the gauge sector, which means wc should consider the 
transformation between the 4 doublets only (i.e., in fia- 
vor space but not in gauge space), including the mix- 
ing between and i/i^. In Majorana fermion represen- 
tation, it is obvious that the allowed flavor transforma- 
tions form 50(8) group. The Lorentz transformations 
*7o7i) *7o72, *7i72 generate SO{3) group in Euclidean 
space, and we also know that 5^(4) = SO{5). The flavor 
symmetry 50(5) and Lorentz symmetry 50(3) actually 
commute. This 50(5) is the largest continuous subgroup 
of 50(8) which can commute with 50(3) and has no 
common element with 50(3) except for identity. There- 
fore 5p(4) is the largest continuous global symmetry. 

If we introduce Nf flavors of fermions, it turns out that 
the emergent symmetry group is 5p(4iVy). We should 
also include the Lorentz symmetry. Here by Lorentz 



Group we mean the continuous group 50( 2,1), gener- 
ated by 7172, 7o7i; 7o72- Note that the physical lattice 
rotation is not identical to the rotation element in this 
emergent Lorentz group. For example, according to Ta- 
ble I, the rotation on lattice -Rgo = Pxy o Px is given 
by 

= ^(71 - 72)(73 + 75)71 
= ^(1 +7i72)(73 +75) 

= ef7i72_L(^3+^5) 

= Dirac 90° Rotation • ^(73 + 75) (100) 
v2 

in the continuum limit. We can sec that the physical 
rotation on lattice is a combination of the Dirac rotation 
and an element in the 5p(4): ;^(73 +75)- This element 
actually exchanges 73 and 75. 

We should also include certain discrete symmetries 
such as time-reversal T, spatial reflections P^, Py, Pxy, to- 
tal parity —1 and charge conjugation. But we know that 
charge conjugation is related to the spin rotation, which 
is included in the 5p(4); Pxy is related to Dirac rotation, 
Px and element :^(73 +75) in 5p(4); and —1 is included 

in 5p(4) as well, namely e™^^ . Therefore the full symme- 
try group of the low energy effective theory for the SU{2)- 
lincar state is Sp{4:Nf) x Lorentz Group x T x Px x Py. 
Such a symmetry group is certainly nnich larger than 
the symmetry group of the lattice model. (The effective 
theory for the SU (2)-linear state does contain terms that 
violate the Sp{4:Nf) x Lorentz Group. But all those terms 
are irrelevant and have vanishing effects at low energies.) 

One can classify the fermion bilinears according to 
their transformation rules under 5p(4) and Lorentz 
group. It is convenient to rewrite the 16 bilinears in 

terms of ip, then in terms of ^: 



tptp 


= —iptp 


tp-yoip 


= -■iplo'4' 


'iphslbi^ 


= 'ipi'lsls'^P 


'ipijij2ip 


= V'i7i72V' 


tpi'yi'tp 


= —tpi^iip 


•0^72 '0 


= -'ipi'y2tp 


ip^ojiip 


= '070 71V' 
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V'7o72V' 


V'naV' = 






■0175-0 


^loizi^ = 


-V'707.3'/' 


i^ioihi^ = 


—07075V' 


i'hiisi' = 


— 0i7i73'0 




-'0i7275'0 




-■0i7i75'0 




-^«7273V' 



(101) 



Horc if there is no minus sign, it transforms as singlet 
under spin rotation. If there is a minus sign after the 
equal sign, it means the fermion bilinear has a (73 in spin 
space, which in turn means the fermion bilinear trans- 
forms as triplet under spin rotation. Triplet should have 
3 components, but in our 16 bilincars we only included 
one of them (the one along z-axis). And the other two 
are fermion bilinears of form ■0-0 and tp^tp^ . 

One can express all the fermion bilinears in terms of 
^f. In summary, one can organize them as in Table III. 
Notice that there are 10 conserved currents of 5p(4), so 
they all have zero anomalous dimension. 

In Table III, we enumerate all the fermion bilinears. 
But what do they correspond to in our original spin 
model? For example, let us look at a particular fermion 

bilinear \E'7o73(Eil(Xicr^'. We will show that this term cor- 
responds to the spin triplet bond order: (— )*»Si x 8^+3.. 
Let us write this fermion bilinear interaction in terms 
of the lattice fermion operators in a unit cell tpi {i = 
0, 1, 2, 3), as shown in Fig. 2: 



Hi =*7o73 O 1 <T* 



( -1 o\ 

-10 
1 

V 1 y 



V'2 

(102) 



Now we can write down the Hamiltonian -ffsu(2)-iincar + 
OHi on lattice. For example, the interaction between 
site-0 and site-1, in terms of /-operator, is 



a. 



(103) 



which basically tells us that (/of/iT/ = + ^) and 



fyn) = ix-e. 

On the other hand, we can write down the spin oper- 
ator So X Si in terms of / operators. Focusing on the 



^-component: 
((So X Si)J 

= -2i((,4/iT./iV/oi - ./oV/u/It/ot)) 

= -2i ((/o\/it) (/o\/u)* - (/o\/u) (/o\/it)*) 
= -8x^ ^ (104) 

Similarly one can show that ((S2 x 83)^) = 8%^, and 
our correspondence is established. 

In Table IV we list the correspondence between the 
field theory operators and original spin operators. From 
this table we know that ferromagnetic order, triplet VBS 
order and staggered spin chirality order all have zero 
anomalous dimension and their correlation function all 
scale as ^ even after the inclusion of the SU{2) gauge 
interaction. We also know that the Neel order and VBS 
order have the same anomalous dimension which turns 
out to be non-zero after inclusion of the SU{2) gauge 
interaction. We noticed that the same 5*^(4) emergent 
group for TT-flux state which rotates Neel order into VBS 
order was found[28] when ignoring the SU{2) gauge field 
effect. 

Here we should mention the work done by Hermele 
et.al[29], where the J7(l)-linear spin liquid was discussed 
and the emergent symmetry is SU{A), and similar clas- 
sification of totally 64 fermion bilinears was done. One 
can recover their result from our formalism. From our 
formulation of S'?7(2)-linear phase, the [/(l)-linear phase 
can be regarded as a Higgs phase where SU{2) gauge 
field is broken down to Uil). Let us assume the remain- 
ing L''(l)-gauge symmetry is along T3 direction. The only 
things one should add to recover their result are: first the 
gauge invariant transformations not only include those in 
Eq.(98); we should also include 



73 (Xi T3 (X) 1, 75 ® T3 (g) 1, 27375 ® T3 



(105) 



So totally 15 elements, and they form a SU{4) algebra. 
Secondly the bilinear with a T3 in the gauge sector is also 
gauge invariant, as shown in Table V. 



G. The effect of SU{2) gauge interaction on 
ySf7(2)-linear spin liquid 

We know that in the continuous limit, the full la- 
grangian should be Eq.(40). The question is, will the 
SU{2) gauge interaction change the low energy behavior 
of mean-field theory drastically? The answer is compli- 
cated. There are two main concerns: spontaneous chi- 
ral symmetry breaking (SCSB) and confinement. SCSB 
means that fermion mass is generated dynamically. And 
confinement means that no excitation with gauge charge 
can show up in the spectrum, and gauge interaction is lin- 
early confining. If any of these happens, the low energy 
behavior of the system is changed drastically and we say 
that the mean-field state is unstable under gauge fluc- 
tuation. In this case, the mean-field SU{2)-\mea.r state 
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Dirac Scalar (5 elements) 


*1 (gi 1 (gi o-*, *i73 (gi 1 ® 1$, $175 (g) 1 (gi 1* 


Dirac Scalar (1 element) 


*j7375 (g) 1 (g) 1* 


Dirac Vector (30 elements) 


*7o (g) 1 (g) (7^, $7073 (g) 1 (g <T$, $7075 (g 1 (g o-$, $i7i72 (g 1 (g 1$ 

$171 (g 1 (g (T$, $17173 (g 1 Ig (T$, $17175 Ig 1 Ig 0-$, $7072 Ig 1 Ig 1$ 
$172 (g 1 (g <T$, $17273 (g 1 <g <T$, $17275 ® 1 Ig 0-$, $7071 Ig 1 Ig 1$ 



TABLE III: Under Sp{4) and Lorentz group, all fermion bilinears can be classified into 3 groups. A group of Dirac scalar and 
Sp{4) 5-dimenion representation, a group of Dirac scalar and Sp{4) singlet, and a group with Dirac vectors in it. For a given 
group of bilinears, they are connected by S'p(4) transformation in each row, and connected by Lorentz group in each column. 
Totally there arc 36 bilinears. All elements in a given group have the same scaling dimension. In the group of Dirac vector, wo 
actually have conserved current corresponding to each column, totally 10 conserved currents. Those are the conserved Sp{4) 
currents, and they all have zero anomalous dimension after inclusion of the SU{2) gauge interaction. 



Dirac Scalar 


$1 (g 1 (g (T$ 


Ncel Order: (-)*Si 


$173 Ig 1 (g 1$ 
$175 lg> 1 (g 1$ 


VBS: {-yySi-Si+^ 

VBS: {-Y'^Si-Si+y 


Dirac Scalar 


$i7375 (g 1 Ig 1$ 


Uniform Spin Chirality: Si • (Sj+x x Si+a,+j,) 


Dirac Vector 


$70 (g 1 lg> fT$ 


Ferromagnetic Order: Si 


$7073 (g 1 (g <T$ 

$7075 (g 1 cr$ 


Triplet VBS Order: {-fySi x Si+o, 
Triplet VBS Order: (-)*»Si x Si+y 


$17172 (g 1 Ig 1$ 


Staggered Spin Chirality: (— )*Si • (Si+a, x Si+a;+y) 



TABLE IV: The correspondence between fermion field operators and original spin operators. In the group of Dirac vector, only 
the density components of each current is presented, since the other components correspond to the flow of these densities. From 
this table, we know that ferromagnetic order, triplet VBS order and staggered spin chirality order all have zero anomalous 
dimension and their correlation functions all scale as A- . We also know that the Neel order and VBS order have same anomalous 
dimension. 



does not lead to a stable algebraic spin liquid. (Actually, 
wc do not know the low energy properties of the model 
beyond the mean- field theory.) 

This problem is a famous and difficult problem in 

QCD, since both effects arc non-pcrturbative. And those 
two effects are related: if there is a mass gap generated 
for the fermion, then below the mass gap there is effec- 
tively no fermion to screen the gauge interaction and we 
only have pure gauge field. We know that pure gauge 
interaction is confining. So logically SCSB will induce 
confinement. The other way, whether confinement will 
induce SCSB, is not clear yet. 

Usually it is believed that for a SU{N) gauge theory 

coupled to Nf fiavor of massless fermions, there is a criti- 
cal iV^[30-32]. If Nf is smaller than Nj, the system have 
both SCSB and confinement. However if Nf is larger 
than Np there is a conformal invariant IR fixed point. 

In particular, in 2+1 dimension, QCD always has a 
stronger interaction at low energies. But if Nf is large 
enough, the fermion screening effect is dominant and the 
renormahzation group flow will terminate at an IR fixed 
point gl The low energy behavior is governed by 

that fixed point, which is deconfined and has no SCSB. 
When Nf — *■ 00, ^ 0, we are in the perturbative 
region. 

To have a quantitative study, we need work within 
large Nf limit [33]. Before we study the SU{2) gauge 



RGFlow gl 



FIG. 4: The fixed point of many Dirac fermions coupling to 
SU{2) gauge field 



fluctuation, it is worthwhile to mention the U{\) case [8]. 
The main result there was that in the large Nf expansion, 
the model remains gapless and the spin-spin correlation 
function is a power law with an anomalous dimension 7. 
This anomalous dimension 7 is calculable in large Nf ex- 
pansion, and up to the 1/Nf (leading) order, 7 is found 



Now we look at the SU{2) case in detail. We will study 
the SU (2)-linear state in particidar. The low energy ef- 
fective theory is Eq.(40). Technically there are two ways 
to do large- A^/ limit. The first way is a complete renor- 
mahzation group analysis. To have an controlled calcula- 
tion, one first does an e-expansion, then studies the renor- 
mahzation group flow, finds out the IR fixed point, and 
the scaling dimension of operators at that fixed point, fi- 
nally sets Nf to be large. This way is generally accepted 
and the result is thought to be reliable. However the 
calculation is complicated. 

Here we do the large-A^/ calculation in a different 
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Dirac bcalar (15 elements) 




*1 


gi 1 (g (7$, $173 (g 1 ig) 1$, $i75 (g 1 (5 


5 1* 








*1 (gi 


rs (g 1$, $173 (g rs (g cr$, $175^x3 


lg(T$ 










$47375 (g T3 (g (7$ 






Dirac Scalar (1 element) 


*«7375 (g) 1 (g) 1* 




$70 5 


D 1 (g)<T*, 


*7o73 Ig 1 (g <T*, *7o76 (g 1 (g o-*, 


$17172 Ig 1 C 


gi 1$ 


Dirac Vector (45 elements) 


$171 C 


g) 1 (g) <T$, 


$17173 (g 1 ® CT$, $47175 (g 1 (g (T$, 


$7072 eg 1 s 


5 1$ 




*i72 C 




4'j7273 ® 1 ® o"'!', >]/i7275 (g 1 (g cr^*, 


*7o7i Ig 1 S 


5 1$ 






*7o73 ® r3 ig) l*, *7o75 ®tz® 1*, ^17172 


(grs (gcr$ 








*i7i73 5 


D r3 (g 1*, *i7i75 Ig r3 (g 1*, *7o72 


g r3 Ig (T$ 








*i7273 5 


5 rs ® 1*, *j7275 (g rs (g 1*, *7o7i (g rs (g <t* 




Dirac Vector (3 elements) 




*i7o 


(g r3 ® 1*, *i7i ® T3 ® 1*, *i72 ® T3 


(g 1* 





TABLE V: Under S'(7(4) and Lorcntz group, the total 64 bilinears can be classified into 4 groups: A group of Dirac scalars and 
5'[/(4) adjoint representation with 15 elements, a group of Dirac scalar and 5'i7(4) singlet with one element, a group of Dirac 
vectors and SU (4) 12-dimension representation with 3 x 15 = 45 elements, and group of Dirac vectors and SU (4) singlets with 
3 elements. 



+ 



FIG. 5: Gauge propagator at leading order of 1/iV/ expansion. 



way [33]. Taking large- TV/ limit: 

Nf — > 00, — > const. 



(106) 



The fermion contribution to any physical quantity can 
be expanded in systematically. This is just a way 
to organize the summation of Feynman diagrams. For 
example, the leading order term usually corresponds to 
summation of fermion one- loop diagrams. The IR fixed 
point can be found by cancellation of leading correction 
to scaling. We will discuss this in detail soon. Here we 
want to discuss whether this approach and the first ap- 
proach are equivalent. There is no general proof that 
these two approaches are equivalent, but in [33] quite a 
lot of examples are presented and it was found that these 
two approaches are equivalent, for example, in the case of 
Dirac fermions coupling to U{X) gauge field, the scaling 
dimension and fixed point found in the two approaches 
are consistent. Gracey et al.[34, 35] calculated large- A/"/ 
expansion of anomalous dimensions of many quantities 
in QED3 and QCD3, including fermion mass, which is 
gauge invariant. Then they compared the results with 
results from usual MS renormalization+dimensional reg- 
ularization, and found they are consistent. 

First let us take a look at the gauge interaction. From 
Eq.(106), is of order -j^. In Figure 5, all leading or- 
der contributions to gauge propagator are summed to- 
gether. The double line represents the leading order 
dressed gauge propagator. 

One can calculate the gauge-gauge two point function 



in Fig. 5. The result is: 



- 7V/'IVlTr(rV^)A"(fc25^, - k^k,) (- ^ 



V 64A; 487rA 



(107) 



Here 1 is the identity matrix in Dirac spinor space. In our 
S'?7(2)-linear example, Trl = 4. And t" are usual SU{2) 
Pauli matrices. Note that the term is independent 
of regularization scheme, but —j^—r is dependent. In 



different regularization schemes, the coefficients in front 



of are difi^erent. 



The result here — 



487rA 



is from the 



Pauli- Villar regularization. 

One immediately sees that at low energy (fc ^ 0), the 
term dominates. The low energy two point function 
of gauge field is: 



NfT^rlUl{kH^.-k^k.)-^AZ 

^ Nf 2 a 



(108) 



The IR fixed point is found by cancellation of the lead- 
ing correction to scaling. Here it means that the — ^p' 
term and the — term should cancel:' 



1 NfTrl 1 



gl 



2 487rA 
127rA 



N- 



f 



(109) 



The dimensionless coupling g* /A describes the strength 
of gauge interaction at the energy scale A. We sec that 
dimensionless coupling g* /A — » as the cut-off energy 
scale A — > and Nj — > 00, as we expected in the renor- 
malization group flow diagram Fig. 4. Remember that 
if the number of flavors of fermions is small, confinement 
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will happen, renormalization group flow will go to some 
fixed point oi gll K ~ 1 as A ^ and perturbation theory 
breaks down. The message is that the fermion screening 
eff'ect in our large Nf model drives the gauge interaction 
to weak limit at low energies. 

From Eq.(107), we know that the gauge field remains 
massless, and in deconfined phase. So the S'C/(2)-linear 
state is also stable under gauge fluctuation. We can also 
sec that the scaling dimension of gauge field AisdA — ^ 
at leading order. With 1/iVj correction, we expect to 
have dA = l + 0(^). Now we can understand why the 
fixed point g1 is IR stable. Suppose we are slightly away 
from the fixed point: 



Nf 



=1 



1 



Tr If f ] 

2g2 IJ ^VJ p,v\ 



(110) 



By power counting, the scaling dimension of 5^ is 4 + 
O(j^), so it is irrelevant. Therefore at low energy 5g 

flows to zero and g1 is IR stable. 

In the above discussion, we have shown that the 
S'[/(2)-linear mean-field state is stable under fermion self- 
interaction and gauge fluctuation in large Nf limit. So it 
is a stable phase. Here by "stable" we mean that the low 
energy behavior does not change drastically from mean- 
field result. When Nf = x we arc back to the mean-field 
result, but if Nf large but finite, the low energy behavior 
is changed from the mean-field result, but not drastically. 
Below we will see what this change is. 



VI. SPIN-SPIN CORRELATION FUNCTION IN 
S'?7(2)-LINEAR PHASE 



In this section we investigate the spin-spin correlation 
function. In a frustrated spin liquid, there is always a 
strong trend to antiferromagnetic long range order. To 
describe this trend, we can calculate the staggered spin- 
spin correlation function in fermion ip formalism: 

{{-irs^{x)s^{o)) = ^ ({U{x)Um - iUy 

(111) 

Here tp is actually an 8-component fermion operator, with 
both spinor indices and SU (2) gauge indices. 

One can see that in field theory language, staggered 
spin-spin correlation function is nothing but the mass 
operator correlation function. At zeroth order, we have 
the free fermion Feynman diagram (Fig. 6). The stag- 
gered spin-spin correlation function in momentum space 
at zeroth order is 



1 

64 



^^p-Tr [Gai3{p)Gi3a{q - p)] 



128 



(112) 




FIG. 6: Spin-spin correlation function at zeroth order. Cross 
means spin operator insertion. 




FIG. 7: Contribution to spin-spin correlation function at or- 
der of -4- 



The first order in expansion involves diagrams in 
Fig. 7. Note that the double line represents the dressed 
gauge field propagator: 



Nf k ^ 



A;2 



sab 



(113) 



These diagrams will give the spin-spin correlation fimc- 
tion an anomalous dimension. According to Eq.(112), the 
spin-spin correlation function in momentum space scales 
like ^ |fc|. With 1/A^/ correction, the correlation fmic- 
tion scales like ~ jfel^^^'''. Here 7 is called the anomalous 
dimension. The contributions from the three diagrams 
in Fig. 7 give (See Appendix A): 



7 



16 

'n^Nf 



(114) 



We see that in the SU(2)-linear state, spin-spin cor- 
relation function remains gapless and algebraic, but the 
scaling dimensions of physical operators, for example spin 
operator, are unusual due to the gapless gauge interac- 
tion. Actually there is no quasi-particle in this phase. 



VII. PHASE TRANSITION BETWEEN SPIN 
LIQUIDS 

A. Phase transition between SU {2)-linear phase 
and S'(7(2)-chiral phase 

1. The effective field theory for phase transition 

In this section we study the phase transition between 
S'?7(2)-linear phase (A in phase diagram Fig.l) and 
S'i7(2)-chiral phase (D in Fig.l). 

We have argued that the SU (2)-linear state is a stable 
phase. What about the S't/(2)-chiral state in Eq.(25)? 
One can check that the fermion spectrum is gapped 
at mean-field level. As for gauge field, because time- 
reversal symmetry is broken in the S'J7(2)-chiral state, 
a Chern-Simons term is generated after integrating out 
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FIG. 8: S(7(2)-linear phase 




5'f/(2)-chiral phase 



the fcrmions. And that wiU also give gauge field a mass 
gap. So the S'C/(2)-chiral state is a fully gapped state. 
We know quantum fluctuation cannot kill a gapped sys- 
tem perturbatively, so SU (2)-chiral state is also a stable 
phase. 

Now we have two stable phases. According to phase 
diagram Fig.l, there is a phase transition between the 
two phases around J2 = 0.35. It also seems that the 
energies of the two phases connect smoothly on the same 
diagram, which indicates the phase transition may be 
continuous. We will study this phase transition in this 
section. 

Fig. 8 shows the plot of the ansatzs of the two phases 
on lattice. Our notation here is 

S'C/(2)-linear phase: S'C/(2)-chiral phase: 



:,z-\-x+y 



(115) 



In Fig. 8, we also show the phases one fermion gains after 
hopping around a plaquette. In S'C/(2)-linear case, this 
phase for the square plaquette is tt; while in 5'?7(2)-chiral 
case, this phase for the triangle plaquette is |. Those 
phases indicate the fluxes through the plaquettes. After 
a time-reversal transformation, the flux will change sign. 
For S'[/(2)-linear case, we have — tt flux, and — tt differs 
from TT by 27r, thus equivalent. This indicates the SU{2)- 
linear phase respects time-reversal symmetry. However 
for the S'C/(2)-chiral phase, we have — ^, which is physi- 
cally different from ^. So the S'[/(2)-chiral phase breaks 
time-reversal symmetry and the parity symmetry. Other 
than that, one can show that 5'?7(2)-chiral phase has full 
translation and rotation symmetry: 

Physical symmetry of S'J7(2)-chiral phase 

={r„r^,i?9oo,c} (116) 

The two phases involved in the phase transition, 
S'J7(2)-linear and 5't7(2)-chiral states, have different 
PSGs, because even their physical symmetries are differ- 
ent. The phase transition breaks time-reversal symmetry. 

The low energy physics at mean-fleld level of the sys- 
tems can be derived by taking the lattice model into con- 



tinuous limit: 

SU (2)-linear phase: 

Lmean = V'^^T/xV' (117) 

5J7(2)-chiral phase: 

Lmean = "4^ d f_,-f + Cr-0 [17375] V' (118) 

To save notation, we dropped the tilde above the fermion 
operator ip. 

Here we see that the a boson field is driving the phase 
transition. In 5[/(2)-linear phase, (a) = 0; in SU{2)- 
chiral phase, (a) 7^ 0. Therefore to understand the phase 
transition, we need to know the dynamics of the a field. 
Let us include quantum fluctuations of all flelds, the low 
energy effective theory is: 

L=Y,^, {d, ~ ^ay) 7m V-. + ^Tr [fljl^] 



1=1 



+ cn/' [17375]'/' + Tr^{d,,(jf + V{a) 

Zp 



(119) 



Here to have a controlled calculation, we again intro- 
duced Nf flavors of fermion. The first line is the Dirac 
fermion couphng to SU{2) gauge field. The second fine 
is the coupling between the fermion and a boson, and 
the dynamics of a boson field. The potential V{(j) is not 
known yet. Nevertheless we know in S'f7(2)-linear phase, 
the dynamics of a boson gives (tr) = 0; while in SU{2)- 
chiral phase, (tr) ^ 0, and the time-reversal symmetry 
is broken. This is similar to the usual formalism of the 
phase transition of symmetry breaking, except for the 
fact that we have gauge field involved here. 



2. The correct effective theory from PSG consideration 

Eq. (119) is the effective Lagrangian for both phases. 
What is the symmetry that the lagrangian should re- 
spect? It should respect the full symmetry of the lattice 
model. Before the symmetry breaking the SU {2)-\mear 
phase has a symmetry described by the S'C/(2)-linear 
PSG. Thus, the effective theory for the phase transition 
should respect the symmetry described by the SU{2)- 
linear PSG. 

We have shown that, under the S'[/(2)-linear PSG, the 
fermion bilinears transform in the way described by Table 
II. For example, the transformation P^^'^: 



P. 



PSG 



#[7375]^' -#[7375]'0 



and 



P. 



PSG 



(120) 



(121) 



In fact, (T is the average of i'0[7375]V' ^-nd transforms in 
the same way as i-ipijsjslip under the S'C/(2)-linear PSG 
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FIG. 9: The behavior of potential V{a) before(left) and af- 
ter(right) the phase transition from 5'C/(2)-hnear phase to 
S'C/(2)-chiral phase. 



FIG. 10: The leading order a boson propagator in expan- 




(see Table II), thus the term crV'[^7375]'/' is nivariant un- 
der the full 5'C/(2)-hnear PSG. We also see that the other 
three possible couplings a^pip, aipjo'ipj a-nd a^pjij2ip are 
not invariant under the S'?7(2)-linear PSG and hence are 
not allowed the effective theory. 

Similarly, the potential V{(t) should also respect a 
—a symmetry and take a form 



V{a) = + 



(122) 



up to quartic order. There is no cubic term since it breaks 
the PP^'^. 

Here we see that PSG tells us the correct form of low 
energy effective theory: 



N, 



L=J2^^{d^- z<t') 7^V. + ^Tr [flJl^] 
i=i ^ 

— 1 771^ 

+ <Ji'[il3l5]i^ + l^id^af + — (7^ + Xa^ (123) 

At the mean-field level, we already have the picture for 
this phase transition, as shown in Fig. 9. We see that 
when m? > 0, (tr) — 0, we are in the S'U'(2)-linear phase; 
when m? < 0, (a) ^ 0, we are in the 5t/(2)-chiral phase. 
rn^ = is the phase transition point. 

At this level the phase transition is second-order, (cr) 
changes continuously from zero to nonzero. The next 
question is, will this transition be second-order after in- 
cluding quantum fluctuations? 



3. The T breaking phase transition does not belong to the 
Ising class 

To answer the above question, we need to count the 
number of relevant coupling constants at the phase tran- 
sition fixed point in the renormalization group sense. If 
there is only one relevant coupling m^, that means the 
phase transition is indeed second-order, and = is 
the critical point. 

We can estimate the scaling dimension of a field. In 
tree level, power counting gives us d„ — ^. But in large- 
Nf limit, the fermion dressing changes strongly. The 
leading order a propagator in expansion is shown in 
Fig. 10. the a boson two point function at leading order 




FIG. 11: a boson contribution to spin-spin correlation func- 
tion at critical point at the first order of expansion. 



IS 
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(124) 



which indicates that the scaling dimension d^ — I + 
0{^). Therefore the scaling dimension of [daY and 

cr^ are both 4 + 0{^), which is larger than space-time 
dimension 3 and thus irrelevant. As for the gauge cou- 
pling 5, the argument in the end of section VG is still 
valid. So 5g is also irrelevant. Now we can safely say 

2 o 

that the only relevant coupling is , whose scaling 

dimension is 2 + 0{^). 

We can calculate the scaling behavior at the critical 
point where a boson is also massless. For example, we 
again compute the staggered spin-spin correlation func- 
tion. At the critical point, apart from the contribution 
from massless gauge field in Fig. 7, we have the contribu- 
tion from massless a boson in Fig. 11 as well. 

So at the critical point, the staggered spin-spin corre- 
lation function not only receives an anomalous dimension 
from gauge field 7 = — ^l^^ (Eq.(114)), but also receives 
an anomalous dimension from the gapless a boson 7'. De- 
tailed calculation shows that at order of 7' — * 
(See Appendix A). 

The total anomalous dimension ^totai is the sum of 7 
and 7': 



Itotal =7 + 7 = - 
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37r27V 



/ 



(125) 



Fig. 12 shows the change of scaling dimension of stag- 
gered spin-spin correlation function during the phase 
transition. Note that in terms of symmetry breaking, 
this phase transition is quite normal: it simply breaks the 
time-reversal, which is a Z2-like symmetry. In our usual 
understanding of phase transition, the Landau-Ginzburg 
paradigm, the symmetry determines the universality of 
the phase transition. Therefore our usual understanding 
for this phase transition should be a Z2-like or Ising-like 
phase transition. However, our study shows it is not the 
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FIG. 12: Change of scaling dimension of staggered spin-spin 
correlation function during phase transition. 
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FIG. 13: At finite temperature T, system shows crossover 
behavior between S'C/(2)-linear phase and quantum critical 
region (dashed line). But the S'C/(2)-chiral phase and quan- 
tum critical region are still separated by a phase transition 
since there is a physical symmetry breaking (solid line). 



case. Although the phase transition is characterized by 
the breaking of a Z2-like symmetry, it is obviously not 
Ising-like. For example, an Ising like transition is gapless 
only at the critical point, whereas in our transition the 
SU{2)-\mear phase is also gapless. The scaling behav- 
ior is very different from the Ising-universality. Actually 
with different value of Nf , the scaling exponent can have 
infinite number of values, which indicates infinite number 
of universalities. 



The above discussion is at zero temperature, but in 
experiments, people can only measure the system at fi- 
nite temperature. What people can see in experiments 
actually should be crossover behavior between disordered 
phase and quantum critical region, as shown in Fig. 13. 





FIG. 14: 5'C/(2)-linear phase 
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FIG. 15: The 4 fermi points of ^2-linear state. 



B. Phase transition between 5'C/(2)-Iinear phase 
and ^2-linear phase. 

1. Z2 -linear phase. 

In this section we study the phase transition between 
S'f7(2)-linear phase (A in phase diagram Fig.l) and Z2- 
linear phase (G in Fig.l). The following is the ansatz of 
the Z2-linear state (Fig. 14). 



^i,i—x-\ 



=?/T , 



(126) 



The fermion energy spectrum of Z2-linear state at 
mean-field level is characterized by 4 fermi points as 
shown in Fig. 15. 

The low energy effective theory is the massless fermion 
coupling to Z2 gauge field. At mean-field level where we 
do not include gauge fluctuation yet, after taking contin- 
uous limit, we have 

Lrnean =i>d^J^4> + "0 [r]T^ijlj5 + ^^"^1121^] 1p (127) 

Once again we need to argue that the Z2-linear state is 
a stable phase. Here by stable we still mean that the low 
energy behavior is not changed after including quantum 
fluctuations. Similar to what we did in Section V, we can 
discuss the stability of Z2-linear state. 

Through a PSG analysis we can show that the fermion 
bilinear term is not allowed and the fermions remain 
massless after including quantum fluctuation [12]. There- 
fore the PSG protects the masslessness of fermion in the 
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Z2-linear state. In addition, Zi gauge fluctuation is al- 
ways gapped. For pure gauge theory on lattice in 
2+1 dimension, it can be in deconfined phase or con- 
fined phase [36]. Here we have massless fermion coupling 
to Zi gauge field, fermion dressing effect should drive 
the system even more likely to be deconfined. We as- 
sume that the gauge field is deconfined, then the gapped 
Z2 gauge fluctuation should be irrelevant. Therefore Z^- 
linear state is also a stable phase. 

In Eq.(127), ryr^ and t^t^ describe the fluxes through 
the red and blue triangle plaquettes in Fig. 14. Since the 
two fluxes are not colincar, the SU{2) gauge group breaks 
down to Z2 gauge group. In Eq.(127), one can say that 
the T] field is driving the phase transition. When (77) = 0, 
wc go back to SU (2)-linear phase, and when (77) ^ 0, we 
are in Z2-linear phase. 

We notice that r] field is not gauge invariant. After 

a local SU (2) gauge transformation, the direction of 
and in Eq.(127) will be rotated. Thus we should use a 
vector field n to describe the fluctuations of ry. Further- 
more, since the rj in front of 7175 and the rj in front of 
7273 should be able to fluctuate independently, we should 
have two vector flelds, say ni and ^2, to describe them. 
So the gauge invariant way to write Eq.(127) would be: 

Lmean = V'^/il/^V' + V' [(^1 " t) 17175 + (n2 • t) 17273] 1p 

(128) 

where Hi and n2 transform as the adjoint representation 
of SU{2) gauge group. Now we include dynamics of n 
fields. To have controlled calculation, we introduce Nf 
flavors of fermions, too. The following is the low energy 
effective theory of the system: 

Nf _ ^ 



i=l 

Nf 



+ V** [(^1 • i7i75 + («2 • t) 47273] i)i + V(ni, n2) 

(129) 



i=l 



where is the covariant derivative of SU{2) gauge the- 
ory, and the form of potential F(ni,n2) is unknown yet. 
The phase transition is described by a Higgs mechanism, 
ni and n2 are Higgs bosons. When (ni) = (ri2) = 0, we 
arc in the S'C/(2)-linear phase; when (ni) ^ 0, (712) ^ 
and (ni) _L (7X2), we are in the Z2-linear phase. 



2. The low energy effective theory from PSG consideration 

What is the symmetry that the low energy effective 
Lagrangian Eq.(129) should respect? Again it should be 

the full symmetry described by the PSG of the SU{2)- 
linear state. We simply need to review Table II again to 



see how n flelds transform under PSG. For example: 

J.PSG . ^[7175]^' ^ -V'[7i75]V' (^3Q<) 
■'/'[7273]V' ^"[7273] 



(131) 
(132) 



To have the term ip [{ni ■ f) ijij^ + (ri2 • t) 47273] i^i 
Eq.(129) invariant under PSG, we should have: 



rpPSG . 


^[7175]^' - 


V'[7i75]V' 


y 


V'[7273]V' - 


-V'[7273] 


pPSG . 


'0[7i75]V' - 


■> -7/! [7273] 


xy 


'0[7273]'i/' - 


-^'[7175] 



-^PSG . ni 

X ' -* 

1-12 

iPSG . ^1 



V 



n2 

PSG . ni 
n2 



xy 



-Til 
ri2 

Hi 
-H2 

-H2 
-Hi 



(133) 
(134) 
(135) 



In summary, the following three transformations of H 
should be the symmetry of the potential V{ni, 712): 



-m 



n2 



-712 



ni 7i2 



(136) 



which strongly constrains the form of potential V^(ni, 712). 
The only allowed gauge invariant form of V{ni, 712) up to 
quartic order is: 



y (711, 712) =-jm'^{ni + n\) + a{n\ + n\) + b{n\n^) 



+ c{Hi ■ 712)^ 



(137) 



Terms like Hi ■ H2 and [Hi ■ n2)(ni)^ are forbidden since 
they break PSG. 



3. A phase transition with no breaking symmetry. 

One can show that the two phases involved in the 
phase transition, the S't/(2)-linear phase and the Z2- 
linear phase, have different PSGs. In [12], PSG of the 
5'i7(2)-linear phase was labelled as SU2BnO, whereas 
PSG of Z2-linear phase was labelled as Z2Azzl3. But af- 
ter projection, the physical symmetries of the two phases 
are identical. They both have the full symmetry of trans- 
lation, rotation and time-reversal: 

Physical symmetry of SU (2)-linear and Z2-linear states 

— {^X) Ty, Px, Py, Pxyi C} (I'^S) 

We are investigating a phase transition with no break- 
ing of physical symmetry. Here the introduction of quan- 
tum order, or PSG is inevitable. Otherwise we do not 
know what is changed during the phase transition. 

At mean-field level, we already have the picture for 
the phase transition. With different values of coupling 
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FIG. 16; Various Higgs condensed phases with different values 
of a,h and c. In phase I, = I(^i2)| and (ni) 1. (n2), it 

is the ^2-linear phase. In phase II, (ni) = and (n2) 7^ or 
vice versa. It is a f/(l)-linear phase which breaks translation 
and rotation symmetry. In phase III, = |(7^2)| and 

(ni) II {n2). It is another i7(l)-hnear phase which breaks 
rotation and translation symmetry. 



constant in potential V(ni,fi2), the Higgs bosons ni,n2 
may or may not condense. If they do not condense, we 
are in the S'L/(2)-linear phase. If they condense in such a 
fashion that (ni) _L (712), we are in the Z2-linear phase. 
Detailed study of the potential shows that if > 0, 
Higgs bosons do not condense. If m? < 0, Higgs bosons 
do condense, and the way of condensation is determined 
by the value of parameters a, b and c as shown in Fig. 
16. There are three different Higgs condensed phases, 
labelled by I, II and III. 

Our Z2-linear phase is phase I. On the mean-field level, 
the phase transition from the S'[/(2)-linear phase to the 
Z2-linear phase can be described by changing from 
positive to negative, and = is the phase transition 
point. 

The next question is whether this mean-field picture 
survives after inclusion of quantum fluctuations. If there 
is only one relevant coupling constant m^, our mean-field 
picture remains valid, otherwise it would fail. 

We now estimate the scaling dimension of ft field. 
Again the fermion dressing is strong in the large- A^/ 
limit. Similar to our argument for a boson in section 
VII A 3, we know that the scaling dimension of n field is 
dfi = l + 0{-^). Therefore by power counting, the terms 

{d^nf' and (those a, 6, c terms) are both of scaling di- 
mension 4 -|- 0(-^), which are irrelevant. However since 
the Lagrangian Eq.(129) is not Lorentz invariant, calcu- 
lating the correction to would be complicated. 

We have just concluded that there is only one relevant 
coupling m? , so the phase transition is second-order and 
vn? = is the critical point. Although a, 6, c couplings 
are irrelevant at the critical point, they are important 
to determine which Higgs condensed phase the system 
would end up. Therefore they are dangerous irrelevant 



1st Order 

Z -linear Phase U( 1 )-linear Phase 
2nd Order 2nd Order 



Disordered Phase 

FIG. 17: The critical theories are the same along the whole 
solid line. 
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FIG. 18: The change of scaling exponent of staggered spin- 
spin correlation function during the phase transition from 
S'l7(2)-linear phase to Z2-linear phase. 

couplings. This can be seen from Fig. 17. Although the 
critical theories for the phase transition from the SU (2)- 
linear phase to Higgs condensed phases are the same, the 
system may change into different Higgs condensed phases 
depending the values of couplings a, b, c. Different Higgs 
phases separate from each other by first order transition 
boundaries. 

The transition from the SU (2)-linear state to the Z2- 
linear state is a phase transition without breaking of any 
symmetry. What are the changes in physically mea- 
surable quantities during the phase transition? Let us 
think about the staggered spin-spin correlation function 
again. On both sides of the phase transition, the fermions 
are massless so the spin-spin correlation functions are 
of power law. But the values of power are different. 
As shown in Fig. 18, in the Z2-linear phase, since Z2 
gauge field is gapped, spin-spin correlation does not re- 
ceive anomalous dimension. At the critical point, due to 
the existence of the massless Higgs fields, the correlation 
function will have another scaling exponent. 

4-. Phase transition from spin liquids to ordered phases 

Spin liquids phases can also experience phase transi- 
tions into ordered phases. Our discussion of phase tran- 
sition between SU (2)-linear phase and chiral spin liquid 
is an example where the time-reversal is broken, but no 
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space translation or spin rotation symmetries are broken. 
Here we focus on the phase transition from spin Uquids to 
phases breaking space translation or spin rotation sym- 
metries. For example, it can go into ferromagnetic phase, 
anti-ferromagnetic phase or VBS phase. 

Suppose our starting point is the S'C/(2)-linear phase 
(or TT-flux phase). Table IV is very useful. For example, 
a phase transition from SU (2)-linear phase to Neel phase 
can be easily understood as the opening of a mass gap 
iliip. Then because at energy scale below the mass gap 
there is no fermion, we are left with a pure SU (2) gauge 
field which is confined in 2+1 dimension. This phase 
transition then have two ingredients in it: fermion-chiral 
symmetry breaking and confinement. Because of the con- 
finement effect the spinous arc always bounded together 
and do not appear in physical excitations. 

Similarly one can study the phase transition from 
S'J7(2)-linear phase to VBS order and ferromagnetic or- 
der. The features for these phase transitions are all sim- 
ilar: opening of a mass gap and the confinement. 

VIII. CONCLUSION 

In this article, we studied the stability of various spin 
liquids, including gapless spin liquid (also called algebraic 
spin liquid). Spin liquids are defined as the disordered 
phases of a spin system. They have the full space-time 
translation and spin rotation symmetries. Different spin 
liquids may have the same physical symmetry. To un- 
derstand the physics of these phases, first one needs to 
understand why they are different despite they have the 
same symmetry. That means we have to classify quan- 
tum states in greater detail than those achieved through 
usual symmetry group. This is the main motivation to 
introduce the idea of quantum order. The PSG is just one 
attempt to characterize quantum order mathematically. 

We find that PSG is very important in understanding 
the stability of algebraic spin liquid. Without PSG, it 
is hard to understand why fermions can remain mass- 
less after including fluctuations around the mean-field 
state. After considering certain PSG transformations 
originated from lattice symmetry, we find that such PSG 
transformations turns into chiral symmetry in continu- 
ous limit, which guarantees the masslessness of fermions. 
Using this idea, we show the existence of an algebraic 
spin liquid - S'?7(2)-linear state whose low energy effec- 
tive theory is a QCD3 with SU{2) gauge group. The 
spin-spin correlation function is also calculated. We 
find that the S'[/(2)-linear state has a large emergent 
symmetry - 5*^(4) x Lorentz Group x T x x Py (or 
Sp{4N) X Lorentz Group x T x P^ x Py for the large A'^ 
model) - at low energies. The lattice model does contain 
terms that violate the 5*^(4) x Lorentz Group symmetry. 
But all those terms are shown to be irrelevant with the 
help of the PSG analysis. 

We also discussed the continuous quantum phase tran- 
sition between spin liquids. Again, the PSG plays a 



key role here. The first transition we studied is a quan- 
tum phase transition that breaks time reversal and parity 
symmetries, which is the transition between the SU{2)- 
linear state and the S'C/(2)-chiral state. Such a Z2 sym- 
metry breaking transition has a well defined order param- 
eter. However, as one can see from the calculated crit- 
ical exponents, the critical point at the transition does 
not belong to the Ising universality class. It is interest- 
ing to see that even some symmetry breaking continuous 
phase transitions are beyond Landau's symmetry break- 
ing paradigm in the sense that critical properties are dif- 
ferent for those obtained from Ginzburg-Landau theory. 

The second transition that we studied is a continu- 
ous quantum phase transition between the 5J7(2)-linear 
state and the Z2-linear state. The two states have the 
same symmetry. Hence we show that continuous phase 
transitions even exist between two phases with the same 
symmetry [19 23]. 

The third transition that we studied is a continuous 
quantum phase transition between the 5't/(2)-linear state 
and a U (l)-linear state. Such a transition breaks the lat- 
tice translation and lattice rotation symmetry. Amaz- 
ingly, we found that the third transition and the; s(x;ond 
transition are described by the same critical theory with 
the same set of critical exponents. So it is possible for 
transitions between very different states to have the same 
critical point. 

All above phenonmena are beyond Landau-Ginzburg 

paradigm for phase and phase transition. Those discov- 
eries suggest that we need rethink Landau-Ginzburg sym- 
metry breaking approach to phase and phase transition. 
We know that stable phases and critical points can all be 
viewed as fixed point in the renormalization group sense. 
If a fixed point has no relevant operator that is allowed by 
the symmetry (or PSG), the fixed point will represent a 
stable phase. If a fixed point has only one relevant oper- 
ator that is allowed by the symmetry (or PSG), the fixed 
point will represent a critical point between two phases. 
In this paper, we found that one cannot use symmetry 
to characterize all the possible fixed points. New kind 
of order beyond the symmetry description exists. We 
showed that how to use the PSG analysis to capture the 
new physics beyond Landau-Ginzburg symmetry break- 
ing paradigm. 

The phase transitions studied here are characterized 
by a change of quantum order in addition to a possi- 
ble change of symmetry. This is why those phases and 
phase transitions are beyond Landau-Ginzburg paradigm 
of breaking symmetry. We emphasize that quantum or- 
der, or PSG, is necessary to understand the correct low 
energy effective theory and the critical phenomena. Also, 
an experimental discovery of a new critical point (with 
unusual critical exponents) implies a discovery of new 
quantum orders. Thus it is very important to measure 
critical exponents even for seemingly ordinary symmetry 
breaking transitions. 

This research is supported by NSF grant No. DMR- 
0433632 and ARO grant No. W911NF-05-1-0474. 
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FIG. 20: Gauge dressed fermion propagator at first order of 



FIG. 19: The dressed gauge propagator in the leading order 
of large Nf limit, which is nothing but the one-loop correction 
in polarization. 



APPENDIX A: DETAILED CALCULATION OF 
ANOMALOUS DIMENSION 7 



The scaling dimension of staggered spin-spin corre- 
lation function ((— )^S(a;)(S')(0)) is calculated by the 
large- A'^ f expansion of quantum field theory. In our 
formalism, one can show that the staggered spin-spin 
correlation function is just the fermion mass operator 
{ijjtjj{x)ipip{0)) correlation function in the effective the- 
ory Eq.(119). By power counting, the scaling behav- 
ior should be (V'i/'(a;)V''!/'(0)) = but quantum fluc- 
tuations change it into {'ijjtjj{x)tljii'{0)) = —t+t^—, where 
7^^ is called the anomalous dimension of fermion mass 
operator. It turns out that the easiest way of calcu- 
lating 7^^, is not to calculate {■ijjijj(x)ijjtjj{0)) directly, 
but to calculate the correlation function of fermion field 
tp: {'ip{x)'ip{0)), and the three-point correlation function 

(ViV(a=)Vi(y)V'(o)). 

Let us firstly calculate the staggered spin-spin correla- 
tion function in SU (2)-linear phase, where the low energy 
effective theory is Eq.(40). We need to understand the 
gauge interaction. In the large- Af/ limit, the gauge field is 
strongly screened by fermions, and under renormalization 
group the coupling g will flow to an IR stable conformal 
invariant fixed point gt ^ Here A is the UV cut-off 

of our theory. To the leading order of -j^j, the dressed 
gauge propagator is as Fig. 19. Let us work within Eu- 
clidean space and Landau gauge, where the bare gauge 
propagator is: 



gauge propagator can be calculated as 



nab (u\ 



(A3) 



fc2(i + n) V""" 

where if we do Pauli-Villar regularization. 
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At the fixed point, where gi — j^yj^ 
propagator is: 



, the dressed gauge 



nab (k) = 



86' 



ah 



NfTFk 



kfxku 



(A6) 



we now study the fermion correlation function with 
first order correction in expansion, as shown in Fig. 20. 
The dressed fermion propagator is: 



Sij [k) 
where 



fc2 



(A7) 



dq^ + Cf8 



8Cf 



{k + 

log{^^ 



NfTFq 
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Sn'^NfTF 'A 



(A8) 



and 



G;l{x,y) = {A;{x)Atiy)} 



I 



dk^ 
(27r)3 



Jk-{x-y)i 



The bare fermion propagator is: 



(Al) 



fc2 



h h 



{i)i{x)'^j{y)) = 



(27r) 



(A2) 



where i,j label the gauge components. The dressed 



Cf = 



(A9) 



Thus we know that the anomalous dimension of is: 



7^^ 



1 8Cf 



2 Sw^NfTF 



(AlO) 



Then we look at the dressed three-point correlation 
fmiction {■ijjijj{x)ip{y)ip{0)) at the order of as shown 

in Fig. 21. Suppose we fix the momentum of ipip to be 2k, 
while ijj and ijj each carry momentum fc, then the tree level 
three point correlation function will be Gz{2k,k,k) = 
'^2 = p- . From the contributions of diagrams in 
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FIG. 21: Gauge dressed three point correlation function at 
order of -j^-. 



Now we can compute 7^^: 



7^^ =A + B + C - 27V, = C + 27^^ 
32Cf 16 



(A13) 



-o-+--a-o- 

+ -0-0-0- 



We can also calculate the spin-spin correlation func- 
tion at the critical point in a similar fashion. The only 



FIG. 22: The contribution of cr-boson to fermion propagator 
at order of where the double dashed line is the dressed 
cr-boson propagator at leading order. 



Fig.21, we will have the dressed three point correlation 

function: 

G3(2fc, A;, fc) = 1 (^1 + + S + C)log{^ (All) 

where A, B, C are the contributions from each corre- 
sponding diagram. Actually wc know that A + B + C = 
7^^ + 27^. So by calculating A + B + C, we will know 
the anomalous dimension of fermion mass operator 7^^ . 

It is easy to sec that A, B arc just from the dressed 
fermion propagator, which has been calculated above: 
A = B = 27^. New calculation needs to be done for 
vertex correction in C. 



Clog{ 
8Cf 



1' 
/ 



NfTp J (27r)3 (g + fc)2(g-fc)2g 



_8CV_ 



(A12) 






FIG. 23: Contributions of cr-boson to three point correlation 
function at order of 

difference is that the a boson becomes massless at crit- 
ical point and contributes to the anomalous dimension 
of correlation functions. The contribution of a boson to 
fermion propagator and three-point correlation function 
are shown in Fig. 22 and Fig. 23. After similar calculation, 
we find that at the critical point. 
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+ 



(A14) 



where the second term comes from the contribution of 
massless cr-boson. The change of scaling behavior during 
this phase transition is shown in Fig. 12. 
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